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Spin Seebeck effect (SSE) holds promise for new spintronic devices with low-energy consumption. The un¬ 
derlying physics, essential for a further progress, is yet to be fully clarified. This study of the time resolved 
longitudinal SSE in the magnetic insulator yttrium iron garnet (YIG) concludes that a substantial contri¬ 
bution to the spin current stems from small wave-vector subthermal exchange magnons. Our finding is in 
line with the recent experiment by S. R. Boona and J. P. Heremans, Phys. Rev. B 90 , 064421 (2014). 
Technically, the spin-current dynamics is treated based on the Landau-Lifshitz-Gilbert (LEG) equation also 
including magnons back-action on thermal bath, while the formation of the time dependent thermal gradient 
is described self-consistently via the heat equation coupled to the magnetization dynamics. 


Introduction. The spin counterpart of the Seebeck 
effect, the spin-Seebeck effect (SSE) refers to the emer¬ 
gence of a spin current upon applying a thermal bias. 
Aside fundamental interest, SSE is of relevance for a 
variety of applications, including spin-dependent ther¬ 
moelectric devices^. Since its first observation^, SSE 
has been studied extensively both experimentally^^— and 
theoretically^^— and for a variety of systems. In partic¬ 
ular, SSE in ferromagnetic (EM) insulators^ hints on a 
magnonic origin of the spin current^. Magnonic SSE has 
some advantages with respect to charge-carrier-related 
spin current in that magnon spin current propagates over 
a length scale of up to a millimeter— while conduction- 
based spin current is usually much smaller due to spin- 
dependent scattering. The spectral characteristics of the 
magnons contributing to the magnonic SSE are still un¬ 
der debate. In a recent spatially resolved experiment^ on 
the magnetic insulator yttrium iron garnet (YIG), the 
measured magnon temperature Tm was related to the 
short wavelength exchange part of the magnon spectrum 
uj{k). An important observation is that a non-vanishing 
spin current emerges even for equal magnon and phonon 
temperatures Tm = Tph- Note that the standard narra¬ 
tive attributes the emergence of the magnonic spin cur¬ 
rent to the difference between magnon and phonon tem¬ 
peratures. A theoretical explanation to the observed fact 
was given in terms of the long-wavelength dipolar part 
of the magnon spectrum^^. If long-wavelength dipolar 
magnons are weakly coupled to the phonons their life 
time is larger and the magnon temperature deviates from 
the phonon temperature. Thus, dipolar magnons might 
contribute to SSE. This explanation though comprehen¬ 
sible, does not exclude the contribution of the long- 
wavelength exchange subthermal magnons in the SSE. 
As was shown in a recent theoretical paper— constitu¬ 
tive issue in the formation of the magnonic SSE is not the 
difference between magnon and phonon temperatures but 
rather the nonuniform magnon temperature profile lead¬ 
ing to a nonzero exchange spin torque and the magnon 
accumulation effect that drives the magnonic SSE. In 
the present paper we study the contribution of the long- 
wavelength (small wave-vector k ) exchange subthermal 


magnons to the magnonic SSE. We find that contrary 
to the short-wavelength exchange thermal magnons, the 
long-wavelength exchange subthermal magnons do con¬ 
tribute substantially to the SSE. Our predictions are 
consistent with the recent experimental report by S. R. 
Boona and J. P. Heremans^. Magnon modes at thermal 
energies in YIG were found not responsible for the spin 
Seebeck effect. Subthermal magnons, i.e., those at ener¬ 
gies below about 30 ± 10 [K], were found important for 
the spin transport in YIG at all temperatures. To assess 
the partial contribution to SSE of magnons with differ¬ 
ent frequencies a;(fc) we analyze the time resolved SSE 
adopting a low-pass filter, as done experimentally^. In 
time resolved SSE experiments the laser modulation fre¬ 
quency is the relevant control parameter. If it is smaller 
compared to the cutoff frequency of the filter then the 
filter does not affect the spin current, otherwise the low- 
pass filter cuts the spin current (external cutoff). Any 
cutoff observed for modulation frequencies smaller than 
cutoff frequency of the filter is thus intrinsic. In our sim¬ 
ulation we implemented this approach, which aside from 
mimicking the experimental situation bears some advan¬ 
tages against a discrete Eourier transform to obtain the 
spectral dependence of the spin current, as detailed in 
the supplementary materials^S. 

We note, that while our study is limited to thermally 
induced magnonic transport in EM insulators, it can 
be in principle extended to include contributions from 
thermally activated carriers in metals or semiconductors. 
The method would rely however on further reasonable 
inputs such as the details of the spin torque current and 
the rescaled exchange interaction parameters. 


Theoretical background. In a EM insulator, the 
low-lying excitations are spin waves describable by 
Landau-Lifshitz-Gilbert (LLG) equation (for a compar¬ 
ison an implementation based on the Landau-Lifshitz- 
Miyasaki-Seki scheme has also been performed with full 
details and results being included in the Supplementary 
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Materials^^) 
d 


dt 


M (r, t) = — 7 M (r , t) X 


04 ^ 
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where M{r,t) stands for the magnetization vector, 7 
is the gyromagnetic ratio, Hq is the external mag¬ 
netic field, A is the exchange stiffness, Ms is the 
saturation magnetization and a refers to the Gilbert 
damping constant. The LLG equation in the linear 
limit has the following solution + iMy{f,t) oc 

exp • r-|-exp (—awfct). The spin wave disper¬ 
sion relation Wfe = 7 (^Hq + and the damping 

which depends on the wave vector aoJk^ indicates that 
short wave vector k (long wavelength) exchange magnons 
are less damped. Thus magnons have different relaxation 
times. Assuming that the magnon-phonon scattering is 
the main source of damping, the magnon-phonon relax¬ 
ation (thermalization) time reads 


= 

' mp 




, k = 2mrld, n = 0, ± 1 , • ■ 


( 2 ) 


where d is the length of the FM insulator. Let us as¬ 
sume that the left edge of the FM insulator is heated 
up periodically with the modulation frequency Wmod (see 
Fig.p. A time periodic thermal bias is meant to mimic 
qualitatively the action of laser pulses (in real exper¬ 
iment the laser electromagnetic energy is absorbed by 
the sample only partly, however. More details are found 
in the Supplementary Materials^S) . The inherent ther¬ 
mal loses relevant to the experiment are treated self- 
consistently by adopting an additional source term in 
the heat equation (see Eq. (SM-3) in the Supplemen¬ 
tary materials^^) . Thus, the temperature profile imple¬ 
mented in the LLG equation in our case is calculated 
self-consistently. For unraveling the role of magnons 
with different frequencies we follow the idea of using a 
low-pass filter, i.e., a filter with an extrinsic cutoff fre¬ 
quency LOc which detruncates the spin current if the mod¬ 
ulation frequency of the laser pulses exceeds the extrin¬ 
sic cutoff frequency of the filter Umod > ^c- The ratio 
between ’’true” /q and measured l{uimod) spin current 

reads \J 1 T (^^mod/^c) • If ^mod ^ 

the measured spin current is not altered by the filter. 
The decay in the spin current in this case is ascribed 
to the intrinsic cutoff frequencies, which in turn are re¬ 
lated to the magnon-phonon-relaxation-time (see Eq. @) 
= 2TTfT^p. In this way different internal cutoff fre¬ 
quencies can be observed. 

Model and simulations. We model a ferromagnetic 
insulator via a chain of FM cells arranged along the x 
axis (Fig. [1]). The total energy density of the system of 


FIG. 1. Schematic of the contribution of the spin waves with 
different wave vectors ki to the spin current in SSE. The 
temperature of the left edge of the system is varied in time: 
To{t) = Too S' (oimodt), where S (tUmodt) is a rectangular pulse 
with a modulation frequency Lu^iod and levels 0 and 1. r„ and 
Mn represent the temperature and magnetization in each in¬ 
dividual cell, respectively, and are calculated self-consistently 
via the heat and LLG equations (Eqs. (U) and ([S|l). 


N cells reads: 


2A 

n—1 ^ n—1 

where a is the size of the cell, Mn is the magnetiza- 
tion vector of cubic cell and A is the exchange stiff¬ 
ness. We use Eq. ([3]) to model Yttrium-iron-garnet (YIG) 
which has been employed extensively in SSE experiments. 
The effective magnetic field acting on the cell reads 

= -sS: = ("“+■ + 


A Gaussian-white noise ffnit) contribution to the 
effective magnetic field (with a correlation function 
{r]niit)r]mjit + At)) = SnmSijS{At)) accounts 

for thermal fluctuations. Here, (■ • •) means average over 
different realization of the noise, n and m are cell num¬ 
bers and i and j are the Gartesian components. The time 
and site-dependent temperature T„(t) obeys the follow¬ 
ing heat equation 

d rr 7-1 _ ^ Tn+ljl) ~ 2T„(t) -|- Tn+l(t) 

dt^~pC a2 ’ 

with the initial and the boundary conditions 
Tn{t = 0) = 0- n = 0,---N + 1 

To{t) = TqqS {uJmodt) ,TN+iit)=^- 


K is the phononic thermal conductivity, p is the mass den¬ 
sity, C is the phonon heat capacity. Too is the tempera¬ 
ture applied on the left edge and S' is a series of rectangu¬ 
lar laser pulses with the modulation frequency Wmod (see 
Fig.p. For solving the heat equation we implemented 
a Forward-Time Gentral-Space (FTGS) scheme^!. The 
hierarchy of the relaxation times for phonons Tph = 
2 l 

(ifa) ~ 10 [ns] (the number and size of cells 

N = 50, a = 10[nm] and phonon thermal conductivity 
K = 6 [W.m“^.K“^]) and magnons Tmp = 2 ^ ~ I0^[ns] 
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FIG. 2. Spin current at the middle of a chain of 50 FM cells 
versus time and for different extrinsic cutoff frequencies (o-’c)- 
We choose a = 10[nm], Ho = 0.057[T], Too ~ 10[K], a = 0.1. 
The system is heated up periodically with the modulation 
frequencies of cUmod = 27r x 10”^ '’[GHz] and the spin current 
is statistically averaged over 1000 realization of the noise. For 
the blue curve no cutoff frequency is implemented on the spin 
current (cjc = oo) but for the red, green, orange and black 
curves the cutoff frequencies are Wc = 27r x 10°'*’[GHz], cjc = 
27 r X 10"°®[GHz], cJc = 27 r X 10“^ ° [GHz] and = 27 r x 
10~^'^[GHz], respectively. 


(ferromagnetic resonance frequency and phenomenologi¬ 
cal damping constant wq = 10 [Ghz], a ~ 10“"^) allows 
an adiabatic decoupling which amounts, to a first order, 
to plug the obtained phonon temperature profile directly 
in the LLG equation and study so the magnetization dy¬ 
namics self-consistently. 


The magnetization dynamics is governed by a set of 
coupled LLG equations 




1 -I- a' 


rM„(t) X 


fff(t) + Jf/Ut) X Hf{t) 


( 6 ) 


For the numerical integration of the coupled stochastic 
differential equations we utilize the Heun’s metho d^^i^^ . 
The spin current tensor is calculated using the following 


TABLE I. Magnon-phonon relaxation times and the cor¬ 
responding frequencies according to Eq. ([2)) for N = 50, 
a = 10[nm[, Ho = 0.057[T[ and a = 0.1. 


n 

0 

1 

2 

3 

4 

5 

|fe| = 2Tm/Na [10®m“^[ 

0.00 

0.13 

0.25 

0.38 

0.50 

0.63 

'^mp 

1.000 

0.303 

0.098 

0.046 

0.026 

0.017 

H^p/2^ [GHz] 

1.0 

3.3 

10.2 

21.6 

37.7 

58.4 



FIG. 3. Normalized spin current (—J^Liriod—\ middle 

^-'2,rXl07[Hz] ^ 

of a chain of 50 FM cells versus the modulation frequency for 
different extrinsic cutoff frequencies (oJc) with the parameters 
a = 10[nm], Ho = 0.057[T], Too = 10[K], a = 0.1. The spin 
current is statistically averaged over 1000 noise realizations. 
For cjc < 27 r X 10® [Hz] the cascades follow the extrinsic cutoff 
frequencies which are characteristic of Low-Pass filter. How¬ 
ever, for tUc > 27 r X 10® [Hz] the cascades occur earlier than 
the corresponding extrinsic cutoff which are a sign of inherent 
intrinsic cutoff in the system. The arrows show the magnon- 
phonon frequencies (H^p/27r in TABLE [T| for different wave 
vectors evaluated theoretically (Eq. ([2])) and coinciding with 
the appearance of intrinsic cutoff frequencies (cascades) in the 
curves. 


formula: J“ = 

where a = x,y,z defines the spin components of the 
tensor, while n stands for the cite number, EaP'y is the 
Levi-Civita antisymmetric tensor and (■ ■ •) means aver¬ 
aging over the different realization of the nois o^’^^i^®’^^ . 
In our model, because of the particular geometry of the 
system (ID chain aligned along x axis) the only non-zero 
element of the spin current tensor is . For the out¬ 
put signal we implemented a recursive low-pass filter— 
Iout{t) = -I- — At). Here 

and lout are the spin currents before and after filtering 
procedure and cOc is the extrinsic cutoff frequency of the 
filter (see Supplementary Materiala^^) . 

Results and discussion. Fig. [2] shows the spin See- 
beck current as a function of time for different extrinsic 
cutoff frequencies uJc- To each cutoff frequency a cer¬ 
tain color is attributed. The filter with the extrinsic cut¬ 
off frequency Wc cuts the spin current if the modulation 
frequency of the laser pulses ujmod exceeds the extrinsic 
cutoff oJmod > Wc- Thus, the larger the extrinsic cutoff of 
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the filter Wc, the larger is the spin current. This is what 
we see in Fig. [21 For small extrinsic cutoff frequencies 
(see Fig. 131) Wc < Stt x 10® [Hz] the spin current is detrun¬ 
cated extrinsically at modulation frequencies uimod = 
smaller than the first intrinsic cutoff frequency of the 
system = 2'kIt^^ . Therefore, no inherent cutoff is 
observed in this case. However, for an elevated extrinsic 
cutoff ujc we observe a cascade of the inherent intrinsic 
cutoff at the frequencies < Wc- All these 

intrinsic cutoff frequencies are in the subthermal regime 
of the magnon spectrum k < k^ax ~ 10® [w -ij6^ (TA¬ 
BLE |T|. To be confident while heating up the system the 
thermal magnons are also activated, the corresponding 
dispersion relation to our parameters is shown in Fig. |31 
As can be seen, magnons with a broad range of frequen¬ 
cies, beyond subthermal regime are also activated. The 
subthermal regime of the spectrum is shown with a small 
green framed. 



FIG. 4. Spin-wave dispersion relation: The yellow back¬ 
ground shows the absolute value of the discrete fourier trans¬ 
formation of rrix -f inny based on micromagnetic simulations^^ 
for a chain of 50 FM cells under a linear temperature gradient. 
Ti = 0, Tiv = 10 [K], a = 10[nm], Ho = 0.057[T]. For small 
k, the dispersion relation reduces to cjfc = 7 ^Hq + 

The small green frame shows the subthermal regime of the 
spectrum. 


To ensure that our findings are not an artefact of a 
particular choice of parameters/model we performed the 
calculations for various temperature gradients, different 
applied external magnetic fields, and varying lengths of 
the chain (Fig. [5]). In all these cases the intrinsic cut¬ 
off frequencies follow the corresponding magnon-phonon 
frequencies (Eq. (|2|)). Furthermore, the use of a white 
noise for a swift time-dependent heating might be ques¬ 
tioned. Therefore, we implemented the Landau-Lifshitz- 
Miyasaki-Seki scheme^ to account for the back-action 
of the magnon subsystem to the surrounding phonon 
bath (see Supplementary Materials^^) and arrived basi¬ 
cally at the same conclusion that small-wave vectors k ex¬ 
change subthermal magnons contribute substantially to 


the formation of thermally activated spin current. Hence, 
our finings are expected to be of some generalities for 
magnon-driven SSE and associated devices, for the em¬ 
ployed schemes are quite ubiquitous, had proven to be 
reliable for finite temperature spin dynamics, and the 
employed system parameters are generic. 
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FIG. 5. Normalized spin current ( ^^mod —"j middle of 

^■'27rXl07[Hz] ^ 

a FM chain versus modulation frequency for different extrinsic 
cutoff frequencies (ojc). R is the number of realization of the 
noise over which the spin current is statistically averaged. 
The arrows show the magnon-phonon frequencies (r^mp/27r) 
for different wave vectors calculated theoretically (Eq. @) and 
compared to the appearance of the intrinsic cutoff frequencies 
(cascades) in the curves. 
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